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Abstract 

One of the goals of the COMPASS experiment is the determination of the gluon 
polarisation AG/G, for a deep understanding of the spin structure of the nucleon. 
In DIS the gluon polarisation can be measured via the Photon-Gluon-Fusion (PGF) 
process, identified by open charm production or by selecting high px hadron pairs 
in the final state. The data used for this work were collected by the COMPASS 
experiment during the years 2002-2004, using a 160 GeV naturally polarised positive 
muon beam scattering on a polarised nucleon target. A new preliminary result of the 
gluon polarisation AG/G from high p? hadron pairs in events with Q 2 > 1 (GeV/c) 2 
is presented. In order to extract AG/G, this analysis takes into account the leading 
process contribution together with the PGF and QCD Compton processes. A 
new weighted method based on a neural network approach is used. A preliminary 
AG/G result for events from quasi-real photoproduction (Q 2 < 1 (GeV/c) 2 ) is also 
presented. 



1 Introduction 

The COMPASS experiment is located in the Super Proton Synchrotron (SPS) accelerator 
at CERN. For a more complete description of the experimental apparatus the reader is 
addressed to [I]. In 2007, the COMPASS collaboration estimated the quark contribution 
to the nucleon spin with high precision [2], using a NLO QCD fit with all world data 
available. This contribution confirms that approximately 1/3 of the nucleon spin is carried 
by the quarks, as demonstrated by earlier experiments [3]. 
The nucleon spin can be written as: 

± = ±AE + AG + L (1) 

AE and AG are, respectively, the quark and gluon contributions to the nucleon spin 
and L is the orbital angular momentum contribution coming from from the quarks and 
gluons. 

The aim of this analysis is to estimate the gluon polarisation, AG/G, using the high 
transverse momentum (high px) hadron pairs sample. The analysis is performed in two 
complementary kinematic regions: Q 2 < 1 (GeV/c) 2 (low Q 2 ) and Q 2 > 1 (GeV/c) 2 (high 
Q 2 ) regions. The present work is mainly focused on the analysis for high Q 2 . However, 
the analysis for the low Q 2 region is summarised in sec. [6j 

For completeness, the slides of the presentation can be found in [3]. 



2 Analysis Formalism 



Spin-dependent effects can be measured experimentally using the helicity asymmetry 

Aa _ a™ - a™ 

~~2o~ an + a n (2) 

defined as the ratio of polarised (Act) and unpolarised (a) cross sections. "flT an d T4 refer 
to the parallel and anti-parallel spin helicity configuration of the beam lepton (\) with 
respect to the target nucleon (ff or JJ-). 

According to the factorisation theorem, the (polarised) cross sections can be written as 
the convolution of the (polarised) parton distribution functions, (A)g,, the hard scattering 
partonic cross section, (A)cr, and the fragmentation function Df. 

The gluon polarisation is 
measured directly via the Pho- y f 1 hf 
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of the gluon inside the nucleon. 

To tag this process directly in ^ a , ^ b j 

DIS a high p? hadron pairs 
data sample is used to calculate 
the helicity asymmetry. Two 

other processes compete with Figure 1: The contributing processes: a) DIS LO, b) 

the PGF process in leading or- QCD Compton and c) Photon-Gluon Fusion. 

der QCD approximation, namely the virtual photo-absorption leading order (LO) process 

and the gluon radiation (QCD Compton) process. In Fig. [TJ all contributing processes are 

depicted. 

The helicity asymmetry for the high p? hadron pairs data sample can thus be schemat- 
ically written as: 

AH(x Bj ) = -Rpgf al^^{x G ) + R LO D A\°{x Bj ) + R QCDC a^ c A\°{x c ) (3) 

The Ri (the index i refers to the different processes) are the fractions of each process, 
represents the partonic cross section asymmetries, Aa l /a l , (also known as analysing 
power). D is the depolarisation factor 0. The virtual photon asymmetry A\° is defined 

1 Ei e ■ qi 

To extract AG/ G from eq. ([3]) the contribution from the physical background processes 
LO and QCD Compton needs to be estimated. This is done using Monte Carlo (MC) 
simulation to calculate Ri fractions and af^. The virtual photon asymmetry A\° is 
estimated using a parametrisation based on the A\ asymmetry of the inclusive data [5]. 
Therefore a similar equation to ([3]) can be written to express the inclusive asymmetry of 
a data sample, A l {£ 1 . 

Using eq. ()3]) for the high p? hadron pairs sample and a similar eq. for the inclusive 
sample the following expression is obtained: 



-"^The depolarisation factor is the fraction of the muon beam polarisation transferred to the virtual 
photon. 
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The term 7 4 corr comprises the correction due to the other two processes, namely the 
LO and the QCD Compton processes. a%, a 2 , 0i, P2, %c, x 'c an d estimated using 

high pt and inclusive MC samples. 



3 Data Selection 

Data from 2002 to 2004 years is used. The selected events have an interation vertex 
containing an incoming muon beam and a scattered muon. As mentioned in sec. [2] the 
data samples are divided into two data sets: the high pt hadron pairs and the inclusive 
data samples. 

Both data sets have the Q 2 > 1 (GeV/c) 2 kinematic cut applied. Another cut is 
applied on the fraction of energy taken by the virtual photon, y: 0.1 < y < 0.9. These 
cuts described previously are used to select the inclusive sample. 

In the high p? hadrons data sample, besides the inclusive selection, events with (at 
least) two outgoing high px hadrons are considered. These so-called hadron candidates 
must fulfill the following requirement: the two hadrons with the highest transverse mo- 
mentum must have pt > 0.7 GeV/c. This requirement constitutes the high p? cut. All 
these cuts additionally correspond to the high p T sample selection. 



4 Monte Carlo simulation 

Important information to be used in the AG/G extraction is obtained from MC simula- 
tion. In this analysis it is fundamental that the simulation describes well the experimental 
data. Two MC samples were produced: one for the high pt sample and another for the 
inclusive sample, to estimate the terms in the set of eq. (j5J). 
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COMPASS 2004 
High-p , Q z >1 (GeVlcf *MC 




The MC produc- 
tion comprises three 
steps: first the events 
are generated, then 
the particles pass 
through a simulated 
spectrometer using 
a program based on 
GEANT3 © and fi- 
nally the events are 
reconstructed using 
the same procedure 
applied to real data. 

For the first step 
the LEPTO 6.5 \7\ 
event generator is 
used together with a 
leading order para- 
metrisation of the 
unpolarised parton 
distributions. The 
MRST04LO set of 
parton distributions 
is used in a fixed- 
flavour scheme gen- 
eration. This set of 
parton distributions 
has a good descrip- 
tion of F2 in the 

COMPASS kinematic Figure 2: Comparison between data and MC simulations - The 
region. distributions and ratios of Data/MC for: inclusive variables: 

NLO corrections x q2 ^ y ^ g ^. row )_ p or hadrons px (2nd raw). For hadron mo- 
are simulated par- menta (3^ row ^ anc j a i so ^he comparison of MC with LEPTO 
tially by including defau i t tunirig . 
gluon radiation in 

the initial and final states (parton shower - PS). 

The fragmentation is based on the Lund string model [5] implemented in JETSET 
[TO] . In this model the probability that a fraction z of the available energy will be 
carried by a newly created hadron is expressed by the Lund symmetric function f{z) = 
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+ pj_, where m is the hadron mass. 



To improve the agreement between MC and data, the parameters (a, b) in the frag- 
mentation function are modified from their default values (0.3,0.58) to (0.6, 0.1). 

The transverse momentum of the hadrons, px, at the fragmentation level is given by 
the sum of the pt of each hadron quarks. Then the pt of the newly created hadrons is 
described by three steering paramrters JETSET parameters: PARJ(21), PARJ(23) and 
PARJ(24). The default values of these three parameters are (0.36, 0.01, 2.0), and were 
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modified to (0.30, 0.02, 3.5). 

The remarkable agreement of the MC simulation with the data is illustrated in Fig. 
121 this figure shows the data-MC comparison of the kinematic variables: xbj, y and Q 2 
(1st row), the hadronic variables, pt for the leading and sub-leading hadrons, together 
with the sum of p^, i.e. YIPti + Pt2 (2nd row), and the momentum p of those hadrons 
(3rd row), also two comparisons of the Yip? variable one using the COMPASS tuning 
and another using the default LEPTO tuning. In this example, it is evident that the 
COMPASS tuning describes better our data sample than the LEPTO default one. 

5 The AG/G extraction method 

In the original idea of the high p? analysis, the selection was based on a very tight set 
of cuts to suppress LO and QCD Compton. This situation results in a dramatic loss of 
statistics. A new approach was found, in which a loose set of cuts applied, combined with 
the use of a neural network [TT] to assign a probability to each event to be originated from 
each of the three processes. The main goal of this method is to enhance the PGF process 
in the events sample, which accounts for the gluon contribution to the nucleon spin. 

The neural network is trained using MC samples. In this way the neural network is 
able to learn about the three processes in order to be disentangled. A parametrisation of 
the variables Ri, x l and a l LL for each process type are estimated by the neural network 
using as input the kinematic variables: %Bj and Q 2 , and the hadronic variables: pri, Pt2, 
p L1 , and p L2 . 

As the fractions of the three 
processes sum up to unity, we 
need two variables to parame- 
terise them: 0\ and 02. The re- 
lations between the two neural 
network outputs 0\ and 02 and 
the fraction are -Rpgf = 1 — 0i — 

I/a/3-02, -RQCDC = Oi-I/a/3-02 

and -Rlo = • o 2 . 

A statistical weight is con- 
structed for each event based on 
these probabilities. In this way 
we do not need to remove events 
that most likely do not came 
from PGF processe, because the weight will naturally reduce their contribution in the 
gluon polarisation, thus enhancing the sample of events thathave a PGF likelihood. 

The resulting neural network outputs for the fractions are presented in Fig. [3] in a 
2-dimensional plot. The triangle limits the region where all fractions are positive. For the 
inclusive sample the average value of 02 is quite large, which means that the LO process is 
the dominant one. The situation is different for the high px sample, in which the average 
outputs are (01) ~ 0.5 and (02) ~ 0.35. Note also that the spread along 02 is larger than 
along o\. 

This means that the neural network is able to select a region where the contribution 
of PGF and QCDC is significant compared to LO, although it can not easily distinguish 
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Figure 3: 2-d output of neural network for estimation 
that the given event is PGF, QCDC or LO; (left) for the 
inclusive sample and (right) for high pt sample. 
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between the PGF and QCDC processes themselves. 



6 High pt hadron pair analysis for low Q 2 region 
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The reason for splitting the Q 2 range in 
two complementary regions is that for the 
low Q 2 region the resolved photon contri- 
butions are considerably higher (« 50 %) 
than in the high Q 2 region, which contains 
practically only the three processes pre- 
viously mentioned. This means that the 
QCD hard scale is also different for both 
Q 2 regions: for low Q 2 the scale is given by 
the high p? hadrons, while for the high Q 2 
is given by the Q 2 value itself. 

A more complicated description of the 
physics than pure QCD in lowest order 
needs to be included in the MC simulation 
for this case. Therefore the event gener- 
ator used in this analysis is PYTHIA 6.2 
[T2] which covers the physiscal processes 
for quasi-real photoproduction. 

In this analysis the selection is essen- 
tially the same as in high Q 2 region plus 
a slightly strict set of cuts: Xp > 0.1, 
z > 0.1, and 5>t > 2.5 (GeV/c) 2 . The 
data sample in this region is 90 % of the 

whole data for all Q 2 range. The weighting method used in the high Q 2 analysis is not 
applied in this case. 

The MC simulated and real data samples of high p? events are compared in Fig. H] for 
Q 2 , y (1st row), and for the total and transverse momenta of the high p? hadron (2nd 
row), showing a good agreement. 

The gluon polarisation in the low Q 2 region is extracted using averaged values as 
shown by this expression: 



p [GeV/c] 



P T [GeV/c] 



Figure 4: Comparison between data and MC 
simulations - The distributions and ratios 
of Data/MC for: kinematic variables: Q 2 , y 
(1st row). Total and transverse momentum 
p and pt of the high p? hadron (2nd raw). 



(6) 

Rfp is the fraction of events in the high px sample for which a parton / from the 
nucleon interacts with a parton f 1 from a resolved photon; A\ is the virtual photon 
deuteron asymmetry measured in an inclusive sample; A/// (A/ 7 // 7 ) is the polarisation 
of quarks or gluons in the deuteron (photon). 

This analysis was performed using a data sample from the years 2002 to 2004. For 
more details about this analysis the reader is invited to look into ref. [13] . 
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7 Results 



The preliminary measurements of the gluon polarisation in low and high Q 2 regions, using 
data from the years 2002 to 2004, are: 



(AG/G) Iow 
(AG/G) high 



0.02 ± 0.06 (stat . } ± 0.06 {syst) 



0.08 ±0.10 



(stat. 



±0.05 



(syst.) 



with xq 
with xg 



Cong 
o.08±8;g| 



The average of the hard scale, /x 2 , for low and high Q 2 is about 3 (GeV/c) 2 . ig 
is the momentum fraction carried by the probed gluons obtained from the MC parton 
kinematics. The result of the measurement for low Q 2 using data from 2002 and 2003 can 
be found in 



Fig. [5] shows these new values of AG/G 
together with the preliminary value from 
the open charm analysis. Also the figure 
shows the measurements from SMC collab- 
oration, from the high analysis for the 
Q 2 > 1 (GeV/c) 2 region [TJ] and also the 
measurements from HERMES collabora- 
tion, for single hadrons and high p? hadron 
pairs analyses [12]. The curves in the fig- 
ure are the parametrisation of AG/G(x) 
using a NLO QCD analysis done by COM- 
PASS [2] in the MS scheme with a renor- 
malisation scale (/i 2 ) = 3 (GeV/c) 2 . The 
dashed line curve is the QCD fit assuming 
that AG > 0, the dotted line is the QCD 
fit assuming AG < 0. It is seen that both 
results from high p? analyses, for high and 
low Q 2 regions, are compatible with each 
other and also, within their xg region, in 
agreement with the NLO QCD fits. 

8 Conclusions 



COMPASS, open charm, prel., 02-06 
COMPASS, high p , Q 2 <1 {GeV/c} 2 , prel., 02-04 
COMPASS, high p , 2 >1 {GeV/c} 2 , prel., 02-04 
HERMES, high p^, all G 2 
HERMES, single high p hadrons, all Q 2 , prel. 
SMC, high p t , Q 2 >1 (GeV/c) 2 
fit with AG>0,[i 2 =3(GeV/c) 2 
fit with AG<0,[i 2 =3(GeV/c} 2 



O) 1 
D)0.8 

< 0.6 
0.4 
0.2 
-0 
-0.2 
-0.4 
-0.6 
-0.8 
-1 



Figure 5: Comparison of AG/G measure- 
ments from COMPASS [2], SMC [Hj, and 
HERMES [15] . The two curves correspond 
to the parametrisation from the NLO QCD 
analysis in the MS scheme with scale at 3 
(GeV/c) 2 . 




The preliminary values of the gluon polarisation for low and high Q 2 regions were pre- 
sented. The gluons were probed at an average scale (/i 2 ) ~ 3 (GeV/c) 2 . Both measure- 
ments show that the gluon contribution to the nucleon spin for xq ~ 0.1 is compatible 
with zero. In that region of xq the presented measurements are in agreement with other 
well known results. 
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